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Abstract: This work sets up a general theoretical framework to study stability of models 
with a warped extra dimension where scalar fields couple minimally to gravity. Our 
analysis encompasses Randall-Sundrum models with branes and bulk scalars, and general 
domain-wall models. We derive the Schrodinger equation governing the spin-0 spectrum 
of perturbations of such a system. This result is specialized to potentials generated using 
fake supergravity, and we show that models without branes are free of tachyonic modes. 
Turning to the existence of zero modes, we prove a criterion which relates the number of 
normalizable zero modes to the parities of the scalar fields. Constructions with definite 
parity and only odd scalars are shown to be free of zero modes and are hence perturbatively 
stable. We give two explicit examples of domain-wall models with a soft wall, one which 
admits a zero mode and one which does not. The latter is an example of a model that 
stabilizes a compact extra dimension using only bulk scalars and does not require dynamical 
branes. 
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1. Introduction 

Extra dimensions are a plausible extension to the Standard Model (SM), and a vast ar- 
ray of specific realizations of models with extra dimensions exist. String theory provides 
a lot of the motivation to include extra dimensions, and there are many string-inspired 
phenomenological models, for example Refs. |l|, ^, ^, as well as purely field-theoretic 
constructions, for example Refs. 0, §• 

Much work has been devoted to the study of the type I Randall-Sundrum (RS) 
model Q, where the electroweak hierarchy can be naturally generated by the warping 
of spacetime in the extra dimension. Such warping is a generic feature of models with 
an extra dimension: the warp factor of the background metric is sensitive to changes in 
the energy density as one moves along the extra dimension, and so non-trivial background 
sources will generically induce a warped metric. In the simplest 5D set-up, the line element 
of a warped space is 

ds^ = e-'^''^y^r]^^dx^'dx'' + dy^ , (1.1) 
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where y is the extra dimension and cr(y) is the warp factor. Scalar fields with profile </>i(y) 
and fundamental branes with tension Aq, located at Ua will provide a source of energy 
density and act to drive the warp factor: 

a"{y) oc <p'.iyf + E ^-^(y - y-^ ■ (1-2) 

i a 

Here, i labels the scalars and a labels the branes. Note that scalars can only provide 
'positive warping', and only if they have a non-trivial profile in the extra dimension. Branes 
can provide both 'positive warping' and 'negative warping' since Aq, can be of either sign. 

In the original RS set-up, the extra dimension was compactified on a circle with y iden- 
tified with —y. The two branes at the orbifold fixed points were required to have tensions 
of equal magnitude but opposite sign in order for the warp factor a to join up consistently; 
the positive and negative warping needed to balance.^ The type II RS model had the 
extra dimension infinite in size, and required only a single brane of positive tension at the 
origin. The warp factor is driven to ±00 at large distances from the origin, localizing 4D 
gravity to the brane. Since this model contains only a positive tension brane, one can 
replace it by a scalar field with a suitable y-profile. Such domain-wall models also localize 
4D gravity and form the basis of a large area of research, see for example Refs. 0, 

Recently there has been interest in a new type of RS-like warped spacetime: a compact 
spacetime where the negative tension brane is replaced with a physical singularity. These 
soft- wall models were originally designed to yield linear Regge trajectories in the context of 
the AdS/CFT correspondence but have since been the basis of actual models beyond 
the Standard Model (BSM) ^ ||, |l6|, 0, |l8| , and also provide a holographic 



dual description of unparticle models p9| , |20| . Our ultimate aim is take the soft-wall model 

and go one step further by removing the final, positive-tension brane at the origin, and 

replacing it with a suitable scalar field profile. The model will then describe a compact 

extra dimension without the use of fundamental dynamical branes, in other words there are 

no branes with tension and localized potential terms in the Lagrangian. This is interesting 

because it is a purely field theoretical construction and requires no appeal to string theory. 
2 

Such a step is not as straightforward as it may first seem, as it is crucial that one 
ensures the stability of the scalar field configuration. Part of this involves stabilizing the 
size of the compact extra dimension. For the RS model, one can utilize the Goldberger-Wise 
mechanism [21, ^| where a bulk scalar field stabilizes the set-up due to localized potentials 



on each brane. Here, the presence of branes and the ability to have a 4D potential for the 
5D scalar localized to the brane is necessary for stabilization to work. Soft-wall models 



can also be stabilized in a similar way |17]; here, one needs a 4D potential localized to the 



^The model also required a bulk cosmological constant to keep cr'{y) constant between the two branes. 
This constant needed to be fine-tuned against the tension of the branes, a generic feature of warped extra 
dimensions with 4D Poincare slices. 

■^The examples we will consider in the following sections, however, will have orbifold fixed points at the 
origin. These are properties of the geometry and one does not need to appeal to string theory to construct 
them. 
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brane at the origin. In order to build a soft-wall model without a brane, we shall need 
an alternative to the Goldberger-Wise mechanism. It is the aim of the present paper to 
show that, with the correct type of scalar background, one does not need any additional 
stabilization mechanism. 

The techniques developed in this paper have wider application than just soft-wall 
models without a brane. The problem of stability of non-trivial scalar backgrounds is an 
import one, and in the case of a single scalar without gravity there exists a constructive 
way of finding the lowest energy, stable solution |23, 24 1. For the case with gravity, things 



are not so simple, but some attempts have been made 26, 27|. The problem is difficult 
because, as discussed above, non-trivial scalar profiles always induce warping of the metric, 
and, furthermore, modes of the scalars mix with spin-0 degrees of freedom in the metric. In 



the literature to date, all successful analyses have relied on the superpotential approach |28| 



also known as the fake supergravity approach |29]. Fake supergravity requires one to choose 



a scalar potential such that it is invariant under supergravity transformations, even though 
the whole theory itself is not locally supersymmetric. In practice this is very simple to 
achieve: one just needs to generate the full potential from a more primitive superpotential. 
Using this fake supergravity approach, one can easily solve the background equations of 
motion, including Einstein's equations. It has been shown, using such a construction, that 
models with an arbitrary number of scalars are free of tachyonic spin-2 modes |^^. For 
a single scalar field, it has also been shown that this construction leads to models free of 



tachyonic modes in the spin-0 sector [30|. In this paper we extend this work to the case of 
an arbitrary number of background scalar fields. We find that the same result holds: there 
are no tachyonic modes in the spin-0 sector. Restricting ourselves to orbifold spaces that 
are symmetric under parity, we also show that set-ups where the scalar profiles are all odd 
do not contain spin-0 zero modes. 

The paper is organized as follows. In Section ^ we discuss the gravity-free case as a 
warm-up and introduction to the superpotential approach. We show that a scalar field 
configuration has a mode spectrum which is always free of tachyons, and which always 
contains a zero mode of translation. In Section |3| we describe the model with gravity 
and the background solutions in the fake supergravity approach. The spectrum of spin-2 
perturbations in this background are derived in Section ^. The usual result of a massless 
4D graviton is reproduced. Section ^ contains the main results of this paper, and shows 
that, in the fake supergravity approach with an arbitrary number of background scalars, 
there are no tachyonic spin-0 perturbations. This analysis is valid for an extra dimension 
of general topology. In Section ^ we specialize to orbifold constructions with definite parity 
and discuss a criterion for determining the existence of spin-0 zero modes. We then analyze 
two specific models: one which supports such a zero mode, and one which does not. The 
latter is an example of a model with a stabilized compact extra dimension that does not 
require any dynamical branes. We conclude in Section ^. 

2. Gravity- free case 

As discussed in the introduction, non-trivial scalar profiles always induce a warped metric. 
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and so realistic models must include gravity. Nevertheless, the gravity- free case is inter- 
esting to study before moving to the case with gravity, and serves also to introduce our 



notation. For one scalar field the general analysis has been completed in Refs. [23, 24 1 and 
for two scalar fields, using the superpotential approach, in Ref. In this section we 

show that, for N scalar fields without gravity, the superpotential approach yields solutions 
which are globally stable, for given boundary conditions, and up to zero-mode translations 
of the configuration. 
The action is 



/ 



d'^x dy 



2 



(2.1) 



where the sum is over 1,...,A^, and there are scalar fields. The potential V is an 
arbitrary function of these N scalars. The equations of motion are 

d^'dM^i - Viii^i}) = . (2.2) 

Here we place a subscript i on y to denote partial differentiation with respect to the field 
<I>j. This notation is used heavily throughout the paper; in general, for a function X that 
depends on the scalar fields, we write 



d d d 



Xii kii'^i}) = T^T^ ---T^ — X({<^i}) . (2.3) 



We are interested in finding static background solutions for the scalars, configurations 
that depend only on the extra dimension: $j = (j)i{y)- In what follows, 'configuration' 
refers to the set of N scalar fields taking on the static solutions given by 4>i{y). 

We can restrict ourselves to potentials generated by a superpotential W^({<I>i}) such 

that 

vim) = Y.hw^({'^^})]' ■ (2-4) 



The subscript i denotes partial differentiation with respect to as per Eq. (|2.3|) . With 
this particular choice of potential, the Euler-Lagrange equations of motion, Eq. (|2.2| ), are 
satisfied so long as the static configuration (pi{y) solve the first-order differential equation 

(l^'M = W,{m) . (2.5) 

A prime denotes a y derivative, and we are evaluating Wi at 4>i{y). Note that to show 
this solves the equations of motion, and in many places throughout this paper, we use the 
relation 

dy ^ d^i dy ' ^ ' ' 

It turns out that such background configurations are always globally stable, up to 
zero-mode translations. To prove this statement, we will look at three things: 1) the 4D 
energy density of the configuration; 2) perturbative modes; 3) the zero mode. 
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We compute the 4D energy density of the system by integrating the 5D stress energy 
over y. We shall assume that V takes the form given by Eq. ( |2.4| ) and that the scalar fields 
are functions of y only, given by (j)i{y). In particular we do not require the scalar fields to 
satisfy Eq. (|2^). We obtain 



E = I Toody 



dy 

= ldW + ^Y.I [^■-Wr{m)]'dy. (2.7) 

i 

The first term in the last line here is just a surface term; it is difference of W evaluated at the 
boundaries of the extra dimension, and depends only on the values of <j)i at these boundaries. 
Given a particular W, which completely determines V, and choices for the boundary values 
of the total energy density of the configuration is minimized precisely when Eq. ( |2.5| ) 
is satisfied. Hence, for given boundary conditions, the superpotential approach yields 
configurations which globally minimize the energy. 

Since the configuration minimizes the energy, the solution must be perturbatively 
stable, and we can demonstrate this explicitly. This will be useful as a precursor to the 
case with gravity. Expanding in linear perturbations ipi about the background 

Mx'',y) = Uy) + v^ii^'',y), (2-8) 

the equations of motion, Eq. ( |2.2| ), reduce to 

-(^'/ + (WikWk, + W,jkWk\^ ipj = Oipi . (2.9) 

Here, repeated indices are to be summed over, and ^lbg indicates that X is a function 
of $j and is to be evaluated using the background solutions 4>i{y); that is, X{{(j)i}). Per- 
forming a Fourier transform, — )• E(pi, we obtain a set of A^, coupled, time-independent 
Schrodinger equations. We can write these Schrodinger equations as [M| 



{dy5ik + Wik\^^){-dy5kj + Wkj\^^ipj = Eipi , (2.10) 

and so, using the results of Appendix we see the system admits solutions only with 
E > {). This is true so long as the perturbations vanish at the boundaries of the extra 
dimension which is a valid assumption. 

Now let us consider the existence of the zero-mode perturbations with E = Q. For 
such a mode to exist, Eq. ( |2.10| ) implies that tp'- = Wijipj, which admits the solution 
(pi = aWi = a(j)'j^, with a a real, non-zero constant. This is the familiar result that the zero 
mode of translation is the first derivative of the background configuration. Note that when 
there are multiple background fields, the zero mode is the mode where all A^ perturbations 
(fi are excited simultaneously with profiles proportional to (f)[. 

As a quick example, consider W = — -y/A/2(i?i>^/3 — v'^ij)), which yields the familiar 
potential V = (A/4)((/)^ — t;^)^. Solutions are cj) = ztw and the well-known kink: cf) = 
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V tanh 



To get the anti-kink solution, one needs to start with — W instead. 
This demonstrates the fact that the superpotential encodes for both V and the static 
solution, and that different superpotentials can yield the same potential. Also noteworthy 
is the fact that, for this choice of ^ = is a solution of the equation of motion but cannot 
be obtained from any superpotential. This is because the cf) = solution is unstable, and 
we have shown that the superpotential approach always yields globally stable solutions. 

3. Warped Background Configuration 

We now come to the main topic of the paper and consider a general 5D theory with gravity 
coupled minimally to N scalar fields, including the possibility of fundamental brane terms. 
The corresponding action is given by 

S = j d^x dy [M^'R + /:„,atter) - a/=^A] . (3.1) 

We are using a ( — h + + +) signature and M is the 5D Planck mass. The scalar-matter 
Lagrangian and brane terms are given respectively by 

>Cmatter = ~\Y. d'"'' ^M^i^N^i ' V{{^i}) , (3.2) 

i 

A = A({$,}) = ^ Kii^imv - Va) . (3.3) 

a 

The sum over i is from 1 to A'^, and y({$i}) is the potential that in general depends on all 
N scalars. The subscript a indexes the branes, ya are their locations, and Aa({$i}) the 
scalar potentials localized to the branes, which includes the brane's tension. Our aim is 
to study the general stability conditions for non-trivial background configurations of this 
model. In this section we discuss background solutions, and then in the following two 
sections analyze spin-2 and spin-0 perturbations. 

For 4D Poincare slices, the background metric ansatz is 

ds^ = e-'^''^y\^,^dx^'dx'' + dy'^ , (3.4) 

where u index the 4D subspace. Einstein's equations are 

Gmn = 2M^^^^ ' 

Assuming the scalar fields only depend on the extra dimension y, and denoting such 
a background configuration by ^i{y), the Euler-Lagrange and Einstein's equations yield 
N + 2 equations for A" -|- 1 functions: 



1 (3.6) 



12a'' - 3a" = [2V ({</>,}) + A ({</.,})] 

</.f-4a'</)--V-({</),})-A,({0,}) = O 
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The first two equations are in fact the sum and difference of the (//i^) and (55) components 
of Einstein's equations. The notation V{{(j)i}) means that the potential is to be evaluated 
with the background fields (pi. Also, we place a subscript i on V and A to denote partial 
differentiation with respect to the field as per Eq. ( |2.3D . It is easy to show that one 
of these equations can be obtained using two others, therefore the equations are not all 
independent. 

Solutions to the above equations can be obtained using the so-called fake supergravity 
approach. In this approach a superpotential W is introduced so that 

i 

where Wi is obtained using Eq. ( p. 3D . It was shown in Ref. |^ that the system of equations, 
Eqs. ( p.6[ ), can be written in terms of the superpotential W and solutions are given by 

However, what we are interested in here are not the solutions of the system of equations. 
Instead, assuming solutions exist, we would like to see if the system is stable or not. 

When studying the gravity-free case in the previous section, we were able to obtain 
the result that the superpotential approach yields globally stable configurations, Eq (^). 
Unfortunately, for the case with gravity included the corresponding analysis does not tell 
us much. The 4D energy density is (recall that we do not assume the relations in Eq (|3.8|)) 



dy 



E = / Toody 



+ 2 1 e-^'^Wim) 



dy . (3.9) 



The first term here is just a surface term, which generally vanishes for a Randall-Sundrum- 
like configuration. The second term is a non-negative integral which is minimized when 
the scalar fields obey the first order fake supergravity equation (as in the gravity-free case). 
The final term is the difficult term. Although it vanishes when a satisfies its first order 
equation, it is not obvious that this minimizes E. So we cannot conclude that the fake 
supergravity approach yields a globally stable background configuration. 

To proceed we shall study local, perturbative stability of a configuration with gravity 
coupled to N scalars. The spin-2 and spin-0 fiuctuations of the metric and the scalar fields 
are treated separately in the following two sections. For the initial stages, our analysis 
will be for an arbitrary scalar potential V. Later on we will need to specialize to the fake 
supergravity approach. 
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4. Spin-2 Perturbations 



The general ansatz which takes into account both spin-0 and spin-2 perturbations is 
ds^ = e-2-(f) [(1 - 2F(x^ y)r/^, + v)] dx^dx'' + [1 + G(x^ y)]^ dy^ , 

Here, we have chosen to work in the axial gauge, /i^s = 0, with transverse traceless part 
d^h^y = r/^'^/i^jy = 0. The Gij (spatial) components of Einstein's equations enforce G = 2F 
which we take from now on. The Einstein tensor has the following non-zero components 

G^^ = 3e'^^T]f,^ {2a'^ - a" - F" + 6a' F' - Ua'^F + 6a" F) 

+ (- le^'^DV - l^M. + ^^'h'^^u + 3(2^'^ - a") V) 
G^5 = 3d^{F' - 2a' F) , 
G55 = 6a'^ - Se^^'^DF + 12a' F' . 
The stress-energy tensor is 

l(l-6FMf-{V + X\^ + (2y + 4A)|bgF- {V, + Xi\^ - cp'.ip', 
1 



(4.2) 



^i'P'i)'-{V + X\,-(.Vi + Xi\^iPi 



J)i5 = d^{(t)'ifi) , 

= ^(<A0' - ^Ibg - 4 V^lbg^ - y^kgV^^ + <Pi^i ■ 

(4.3) 

Here, repeated i indices are to be summed over. Taking the 4-trace of the (/Ui^) part of 
Einstein's equations shows that the spin-2 and spin-0 perturbations completely decouple 
from one another. See Ref. [^] for a discussion of the gauge degrees of freedom and 
the decoupling of spin-2 and spin-0 sectors. The spin-2 perturbations are the easiest to 
analyze, and we look at them first, returning to the spin-0 sector in the following section. 
The Euler-Lagrange equations for the N scalar fields do not contain h^i,; these equations 
are given in the next section. 

Using the background equations, Eq. (|3.6D , the equation for /i^^^ is 

-e2-nV-/i';^ + 4c7'/i;, = 0. (4.4) 

There is always a zero mode, /i^j^ = 0, which is normalizable. It is the well-known 4D 
massless graviton [^. In conformal coordinates z defined by dy = e~'^ dz with rescaled 
h^v = e^'^/'^h^u we have 



-K. + [ja'^ - -a"j = Uh^, . (4.5) 
This is a Schrodinger-like equation. It can be rewritten in a self-adjoint form as 

(5. + Sh){-d, + Sh)h^, = DV , (4.6) 

with 

Sh = -^a' . (4.7) 
Appealing to Appendix we see that there are no tachyonic modes in the spin-2 sector. 
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5. Spin-0 Perturbations 



In the previous section we have seen that the spin-2 sector decouples from the spin-0 sector, 
independent of the number of scalar fields involved. We now consider spin-0 perturbations 
around the background solutions cr(y) and (piiy). The metric ansatz is a restricted version 
of Eq. Q with /i^^, = p 



-My) 



1 - 2F(x^, y)] rj^^dx^'dx" + [1 + 2F{xf',y)f dy^ 



(5.1) 



^^{x^',y) = (j)i{y) + ^i{x'',y) . 

We work to linear order in the perturbations F and ipi. The equations for these perturba- 
tions consist of two of Einstein's equations and the Euler-Lagrange equations: 

6M^{F' - 2a' F) = <p[ip, , (5.2) 

6M3(-e2-nF - 2a' F' + F") = 20',(^', + 2 Aj^g F + A^g , (5.3) 

e'"^n^, + ^1 -4a'ip',-6F'cP',- m + 2X,\^F- iV^J + = . (5.4) 

Note that one of the redundant Einstein's equations has been omitted and it can be easily 
obtained by using Eq. ( |5.2[) and Eq. (|5.3D . 

We now go to conformal coordinates z to eliminate the e^"^ factor from the D's, and 
rescale the fields by a factor of e^*^/^ to help eliminate first derivatives. That is, we make 
the following change of variables and field redefinitions 

dy = e-'^^y^z , a{y) = a{z{y)) , Uv) = M^v)) , 



F{y) = e2 



F{z) 



<^{^(s/)) 



(5.5) 



In terms of the conformal coordinate z and the new fields F and c^j, Eqs. (5.2) through (| 
become, respectively. 



6M^{-WF + F') 



6M^ 



-□F + F" + 2a' F' + (|a'2 + fa") f] = 20',(^', + 2e~^'' X\y^ F 



+ {3a'd + e-^'' Ail 



bg 



Oifi + (p'l + 



9 /2 , 3 //\ ~ 

-fa + fo- )ipi 



?>4>i (3a'F + 2F' 



e-'" m + 2A.)|b„ F - e-^-^ (y,, + A,,-)lb, = . 



(5.6) 
(5.7) 
(5.8) 



We can make judicious use of Eq. ( |5.6| ), and its derivative, to eliminate all first derivatives 
of F and (pi in the other two equations. We also make use of the background equation for 
6~. This allows us to obtain 



6M3(F' - \a'F) 
6M^ \dF + F" - 



0, 



□(^, + + + ^a") + (-40^' + 2e 



'2u 



+ 2e-2- AlbgF + ( -2(/>^' + e-^" A^ J = 



F 



-2a 



This is a set of + 2 equations for + 1 functions of z. 



. 



(5.9) 
(5.10) 

(5.11) 
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5.1 Simplifying the perturbation equations 

The main aim of this section is to manipulate and simplify the above A'^ + 2 equations 
in order to obtain a set of coupled Schrodinger-like equations, which can be solved to 
determine the spin-0 spectrum. Before simplifying the equations, we need to understand 
precisely their redundancy. Let a, (3 and 7j be equal to the left-hand-sides of Eqs. ( ^^ ) 
through ( p. 11 ) respectively. Then the Euler-Lagrange and Einstein's equations amount to 



a = /3 = 7* = 0. (5.12) 

One can show that the following relation holds in the bulk: 

Da + a" - (|(t'2 + |(t") a - /3' + \a' ^ + (fi'^i = . (5.13) 

Thus, there is a certain amount of redundancy in the + 2 equations. The redundancy 
is quantified precisely by Eq. ( |5.13| ). In particular, if we solve a = /? = then we auto- 
matically have (p'^'^i = and solving for — 1 of the 7j's is enough to solve for the system 
(that is, we have eliminated one of the scalars (pi). Counting the number of free integration 
constants (ICs), we have 2 for a = and /3 = together, and 2(A^ — 1) for — 1 of the 
7j's. This is a total of 2A^ free ICs. 

We can make other choices for eliminating an equation. For example, if we choose to 
solve all 7i = and /3 = 0, then we will have an equation left over for a: 

Ua + a" - (|cj'2 + la") a = . (5.14) 

Solving /3 = 7j = only solves a such that it is a solution to the above differential equation. 
Counting the total number of free ICs: 2 each from all N of the 7j's and 2 from /3 = 0. 
However, 2 of these ICs have to be used in order to pick the a = solution of Eq. ( 5.1^ ) 



and so we end up with a total of 2A^ free ICs for the system, the same as in the alternative 
choice above. 

We can also choose to solve a = and 7^ = with 2A^ + 1 free ICs. Then Eq. ( 5.13| ) 
becomes — /?' + ^o"'/3 = and we need to have 1 constraining IC from demanding that we 
pick out the /3 = solution from this differential equation. Thus there are again only 2A'^ 
free ICs for the system. 

These arguments hold for solutions of the system in the bulk. Without branes, the 
system is completely specified by 2A^ integration constants. If there are branes in the 
set-up, the situation is the same, as the values of fields on the brane are continuous and 
hence the same as the values in the bulk just to the sides of the brane(s). Derivatives of 
the fields are discontinuous on the brane, and hence undefined. Instead, one has freedom 
to choose the field derivatives on one side of the brane or the other, and then uses the 
'jump' conditions to determine the derivative on the opposite side. The jump conditions 
are found by integrating Eqs. ( 5.10| ) and ( |5.11| ) over each brane; see |22| for details. 



We choose to solve the system /3 = 7^ = 0, as it leaves us with the most 'symmetric' 
system of equations. Solutions to this particular set of equations will include all solutions 
of the true system (by that we mean Eq. (^l2D), as well as some additional solutions which 
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do not satisfy all the boundary conditions of Einstein's equations. For us, this is all we shall 
require, as we are going to show that the extended set of solutions does not contain certain 
modes (tachyonic and/or zero modes), and so the full system cannot therefore contain such 
modes. 

If one makes the definitions 

^i{z) = M^l^^iiz) , (5.15) 
F{z) = x{z)/^, (5.16) 

then the equations /3 = 7^ = 0, Eqs. ( |5.1C1| ), and ( |5.11| ), become 

□X + X" - (Voo + ^oo)x - (Vo. + BQi)^i = , 
□V'* + < - {V^J + 13,j)i>j - (Voi + Boi)x = , 



(5.17) 



where 



Voo = l^'^ + \a" , 

Vo. = -^j^^'l , (5.18) 



and the brane terms are 



Note the symmetry of the cross-coupling in Eq. ( |5.17| ). We can use Eq. ( 5.17| ) to solve 
for the physical spin-0 spectrum, which includes the eigenvalues and corresponding extra- 
dimensional profiles. 

A physical mode is defined as having the same x^^ dependence in x and V'i; but possibly 
different z dependence. Separation of variables then proceeds by defining 

x(x^z) = /o(%(x^, 

i^i{x^,z)=h{z)r^{x^^), 

where frn{z), m = 0,1, . . . , N are the extra-dimensional profiles and rj(x^) the 4D mode. 
To find the spectrum of mass states of r], we perform a Fourier transform on the variable 
x^: Or] —7- Erj, where E is the eigenvalue, proportional to the mass squared of the rj mode. 
Eq. (5.17) now becomes a coupled eigenvalue equation of the form 

(5.21) 





_ / -di + Voo + ^00 Vo, + Boj 
\ Voi + Boi -d%j + Vij + Bij 

In the matrix multiplication there is a sum over repeated j indices from 1 to N, and we 
shall use such a notation in subsequent equations. 
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This equation is a system of + 1 coupled Schrodinger equations, with a symmetric 
coupling potential. The symmetry implies that the eigenvalues of the system will be real, so 
long as certain boundary conditions are satisfied. This equation is one of the main results 
of the paper. It allows one to determine the spectrum of physical spin-0 modes with N 
scalars coupled to gravity. It is valid for an arbitrary potential V, with or without branes 
at the edges of the extra dimension, and for generic topology of the extra dimension. 

5.2 The case without branes 

We are now going to specialize to the case where the brane terms are absent; A = 0. A 
particular example of such models is one where a domain wall replaces the fundamental 
brane |^]. In this case we know exactly how the perturbations behave at the boundaries, 
and can proceed to determine the spectrum. We shall show, using the fake supergravity 
approach, that there are no tachyonic modes for these models. We present conditions for 
the non-existence of tachyonic modes for models with fundamental branes in the following 
section. 

In order to prove stability for a system with scalar fields and the gravitational 
perturbation Xi we need to show that the eigenvalues E are strictly positive for the given 
boundary conditions. Our analysis so far has been for a general potential V and general 
background configuration of scalar fields. If we specialize to configurations generated by 
the fake supergravity approach, we can prove that the eigenvalues E of this system are 
non-negative, and in some cases strictly positive. 

The key for proving such positivity is the observation that one can write the pertur- 
bation potential as 

V= h'^A=S' + S', (5.22) 



where 



(5.23) 



Then, 

H = {d, + S^){-d, + S). (5.24) 
We have shown in Appendix ^ that the eigenvalues satisfy E >Q whenever 

#(-a, + S)^|boundary = 0, (5.25) 

where ^ = {foJiV- 

As is obvious from the above equation, boundary conditions play a crucial role on the 
positivity of the eigenvalues and therefore on the stability of the system. Here we discuss 
the boundary conditions with general N for the following specific topologies: 

Full-interval space: On a full interval there are no restrictions, such as parity or 
Dirichlet or Neumann boundary conditions, on the perturbation wave function ^. We 
only require that the perturbations are normalizable over the full space. The boundaries 



- 12 - 



are at zi and Z2, which can be at infinity or at finite values as in the case of a soft- 
wall model, where space ends at a physical singularity. As seen in Eq. ( ^.5D solutions 
in conformal coordinates z are related to the ones in y coordinates through the warp 
factor which diverges at large z. We see that in order to have normalizable solutions, the 
rescaled perturbations in conformal coordinates must vanish at the boundaries. Therefore, 
^"^ {-dz + S) ^'U=2i,2 = is satisfied. 

Half-interval orbifold space: Taking a full-interval space and identifying y with 
—y yields an orbifold, with effective boundaries at z = and z = zi. Note that for such 
a topology there will be a non-dynamical brane at the origin, but this will not contribute 
to the brane terms B in the effective potential for the perturbations. Furthermore, we do 
not need to appeal to string theory for information on the physics of this non-dynamical 
brane; it is just an orbifold fixed plane. The boundary conditions for ^ at z = are 
now equivalent to a choice of parity for the warp factor and the scalars. The warp factor 
must be even in order to localize gravity. Scalar fields that are odd vanish at z = 0, even 
scalars have dz^i vanishing at z = 0. It is also easy to verify that ^'^'S' has odd parity, 
and therefore it also vanishes at z = 0. As such, the condition ^'^ {—dz + S) ^'|2=o = 
is always satisfied. For the z = zi boundary, the scalar perturbations vanish as for the 
full-interval case discussed above due to the normalizable condition on the perturbations. 

Having explicitly shown that the boundary terms vanish, which is expected for a 
Hermitian Hamiltonian corresponding to a physical problem, we can conclude that without 
dynamical branes models with scalars coupled to gravity do not have any tachyonic 
modes. Note that our results so far apply to both of these scenarios while in the following 
sections for explicit examples with two or more scalar fields we will concentrate on the 
half-interval case. 

5.3 The case with fundamental branes 

For completeness we would also like to give sufficient conditions for the non-existence of 
tachyonic modes for models with A^ scalars coupled to gravity in the presence of funda- 
mental branes with A 7^ 0. The profiles for the scalar fields for models with fundamental 
branes satisfy the effective Schrodinger equation given by Eq. ( ^.2lD where the brane terms 
Bmn are given in Eq. ( 5.19 ). Repeating the analysis of Appendix ^ with brane terms we 



find that the effective Schrodinger equation can be written as 

j dz\S^'^ + j dzdz(^^S-^^+ j dz (^^t^^) = E j dz\^^ , (5.26) 



where S = {—dz + S) with S given by Eq. ( 5.23| ) and the brane terms are 



B = E^-"' ^ ' S{z - z„) , (5.27) 

where the sum over a is over all fundamental branes in the model. For £^ > the sufficient 
condition is then modified to be 

*t {-dz + S) ^-Iboundary + > ■ (5.28) 
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Analyzing this requirement for different models with different brane potential terms A" is 
outside the scope of this work. 



6. Analyzing Zero-Modes 

As we have seen in the previous section, the eigenvalues of the spin-0 sector are guaranteed 
to be non-negative when constructing configurations using the fake supergravity approach. 
We must still analyze the existence of zero mode solutions to show complete stability. In 
particular this will ensure that the size of the extra dimension is stabilized. 

In this section, using the formalism we have developed, we will first show that for the 
= 1 case the eigensystem reduces to a single Schrodinger-like equation, and we determine 
the mass gap of the spin-0 spectrum. For configurations with N > 2 scalar fields, analyzing 
the general properties of the mass spectrum is a very difficult problem. As such, we restrict 
our attention to orbifold spaces where the fields have definite parity and discuss a criterion 
for the possible existence, or lack, of zero modes. Guided by this criterion, we analyze 
the zero modes in two example models with specific superpotentials. We shall explicitly 
construct a model with N = 2 scalar fields that does not have a zero mode and therefore 
is stable. 



6.1 The A = 1 case 

In this section we will analyze the system with one scalar coupled to gravity and show 
that it is stable with positive mass eigenvalues. The analysis is valid for general topologies. 
Stability for models with one scalar field and no branes have been analyzed previously in 
Refs. [p^, 3C]. For models with a soft- wall and a fundamental brane stability has been 



proven in Ref. |17] and for RS models with two fundamental branes in Refs. p^, p2|. 



We start by analyzing the Einstein constraint equation, Eq. (|5.9|) , and the second-order 



equation for x from Eq. ( 5.21 ). They are, respectively. 



-/o + 5oo/o + %/i = 0, (6.1) 
-fo+Voofo + Voin = Efo . (6.2) 

As before, repeated i, and later j, indices are to be summed. Combining them in an obvious 
way yields 

-/o + Q c Jo + Voo ^r-^ — ] fo+ \ Vo^ — — fi = Efo . (6.3) 

We can eliminate the /q term, and obtain a Schrodinger-like equation, by defining 



/o = \J SoiSoi g . (6.4) 

This gives 

-g" + (^2 _ ^/ ^ g ^ J-^ BJ, = Eg , (6.5) 



-14- 



where 

_ So iSijSoj 



'-'Ofc'-'Ofc 

Si = Voi 5r-^ . (6.7) 



Equation (6.5) can be used to prove the non-existence of a zero mode for a theory with 



one scalar field. For = 1 we have: 

A = Su, (6.8) 
A = . (6.9) 

Following the arguments in Appendix^, we see that for wave functions ^' 7^ 0, the eigenval- 
ues of this system are non-negative since the boundary terms vanish as we have discussed 
previously. Equation (6.5) simplifies to 



(a, - A){-d, - A)g + So^Soig = Eg . (6.10) 



Now we will use the same trick as in the analysis of Appendix ^ and multiply Eq. ( 6.1C| ) 



from the left by g* and integrate over the extra dimension z. Up to a surface term that 
vanishes we have 



dz 



+ / dz\Soig\'^ = E / dz\g\'^ . (6.11) 



Consider the existence of a zero mode, with E = 0. In order to satisfy Eq. ( 6.1l| ) for 



a generic superpotential W, and for the given boundary conditions by which the surface 
terms vanish, both of the terms on the left-hand side must vanish simultaneously. For a 
generic W this means that the field g, and hence /o, has to vanish. This shows that for 
nontrivial /o there does not exist any zero mode. 

We can push the analysis a little further and provide a lower bound on the mass gap 
to the first spin-0 state. Dropping one of the terms on the left-hand side of Eq. ( |6.11| ) gives 
the following inequality: 

E J dz \gf > J dz\Soig\^ (6.12) 
> min[5^i] / dz\g\^ , (6.13) 



where min[S'Q]^] is the minimum of Sq^ over z. Thus we obtain a bound on the mass of the 
first spin-0 state: 

E > min[5^J . (6.14) 

6.2 The iV > 2 case 

Set-ups with more than one scalar field can in general posses a zero mode. Neverthe- 
less, we shall provide a simple criterion that can be used to find models which do not 
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have a zero mode. Prom now on we must speciahze to half-interval orbifold spaces so that 
the N scalar fields have definite parity. This enables us to formulate the following criterion: 

Zero-mode criterion: For a system of definite parity with N scalar fields that couple to 
gravity, the number of independent normalizable zero modes with E = is at most equal 
to the number of even fields. 



This criterion is essentially a statement about integration constants. In the fake super- 
gravity approach, the background configuration is given by the solutions to the first order 
equations, Eq (|3.8| ). This is a system of linear differential equations for -|- 1 functions, 
which require A^-|- 1 Dirichlet boundary conditions for a unique solution. The value of (t(0) 
can always be chosen to be since any other constant shift can be obtained by a redef- 
inition of the coordinates x'^. The restriction to orbifold models further eliminates those 
integration constants associated with odd-parity scalars, since their field value must also 
vanish at y = 0. Unique solutions to the fake supergravity equations are then parametrized 
by the integration constants of the even-parity fields. The final point is that zero modes 
move us continuously through this space of solutions, so there cannot be more zero modes 
than the number of even fields. 

For our above argument to hold, we must show that the zero modes do indeed take us 
from one solution to the next. Looking at Eq. ( A. 3 ) it is easy to see that a solution with 
E = is satisfied when {—dz + S)^ = 0, that is 

- /o + 5*00/0 + Soifi = , -, ,.x 

, (6.15) 

— /, + Sojfo + Sijfi = , 



where the matrix elements Smn are given explicitly in Eq. ( 5.23 ). Now suppose there are 
normalizable zero mode solutions (/g"'', /f^^) which satisfy the above equations. Working 
in y coordinates, we now add these zero mode solutions as perturbations to the background 
configuration to define the following new fields 



a(y-) = a(y) + 6^/r(y), 



'12 ^° ■ (6.16) 
4>i (y) = 0.(y) + ee3/2-M3/2/i°^(y), 
where the new coordinate y is expressed in terms of the zero mode perturbations as 

dy = dyll + 2e^fS'A . (6.17) 



12 

e in the above equations is a small parameter that parametrizes the small perturbations. 

(0) ' (0) ' 

Solving for the derivatives of the zero mode solutions, /q and f- , in terms of the 
zero mode solutions themselves, by using Eq. ( |6.15D , it is easy to show that the new fields 
defined in Eq. ( |6.16D satisfy 

^-Hy) = 7rLw {Hi}) , 

f (6.18) 
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up to first order in the perturbation parameter e. We see that these new fields are them- 
selves background solutions with the same superpotential W. This demonstrates that the 
zero mode solutions for the system of A'^ scalar fields coupled to gravity translates from one 
background solution to another, taking us from one set of integration constants in the y 
frame to another set of integration constants in the y frame. Note that this result is valid 
up to first order in e, so the zero mode solutions have to be strictly much smaller than 
the background solutions, restricting our reasoning to normalizable zero modes. Since any 
given zero mode is associated with continuously changing a Dirichlet boundary condition, 
the number of physical zero modes cannot be more than the number of even-parity scalars. 
This completes the proof of our criterion. 



6.3 Explicit Examples 

In this section we look at two specific domain wall models with N = 2 scalar fields coupled 
to gravity. We have already shown that the spin-0 spectrum has strictly non-negative mass 
eigenmodes, and proven a criterion relating the existence of zero modes to the parities of the 
background scalar profiles. The two models to be presented will form explicit realizations 
of this criterion: one has an even scalar profile and a zero mode, the other all odd profiles 
and no normalizable zero mode. Aside from this difference, both models are qualitatively 
the same. The field $i will play the role of a dilaton and has a background solution which 
diverges at finite y, generating a physical singularity and cutting off the extra dimension. 



effectively compactifying it. This is known as a soft wall 11, |T^, 17 1. The second 
field ^2 takes the form of a kink, creating a domain-wall whose purpose is to replace the 
positive tension brane in usual soft-wall set-ups. The models we consider are on a half- 
interval orbifold space with definite parity, and the domain-wall sits at the origin, acting 
as an effective boundary of the extra dimension. A domain-wall soft-wall model is an 
appropriate name for this type of set-up. In this section we are concerned primarily in 
the stability of such models, and do not discuss any other phenomenology. Whether or 
not these models solve the hierarchy problem is an interesting question which we intend to 
address in future work. 

Example 1 

The first model we consider has the following superpotential 

iy($i,$2) = (a$2 -6^1) e''*! , (6.19) 

where $i and $2 are the dilaton and the kink fields respectively, a > 0, 6 > and v are 
parameters in the model. We write the fields in terms of their background solutions and 
perturbations as 

$1 (x, y) = (y) +(pi{x,y), 
^2(2;, y) = My) + (p2(.x,y) . 
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Figure 1: Numerical solutions for the backgrounds u, 4>i and 02 for example 1 with the parameter 
values a — 1, b ~ 1 and ly = 1.4. The free initial condition for 0i(O) is chosen to be —0.1. The soft 
wall is at z = 3.41 where cr and 4>i diverge, ending the spacetime at a physical singularity. These 
solutions have definite parity and we have only plotted the z > half. 

The background fields satisfy 

dy 

^^2 = (a-3b4) 6"^^ (6.21) 
dy 

We choose the background solutions such that cpi has even and has odd parity, that is 

^4i\y=o = 0, </.2(0) = 0. (6.22) 

Using our criterion we know that this model will have at most one zero mode solution 
since we have one field whose background solution is even, which is the dilaton field. The 
nice thing about this simple superpotential choice is that the zero mode solution can be 
found analytically even though an analytic solution for the background fields is not possible. 
Fig. 1^ shows the numerical solutions for the background fields for a particular parameter 
space point. Although we are not concerned here with the parameters of the model, an 
important thing to mention is that our choice of = 1.4 allows us to satisfy the equations 
of motion at the singularity; see Ref |17] for details.^ 



Going to the conformal coordinates using Eq. ( |5.5D and using the field redefinitions in 



^Note that, in addition to $i, the parameters of our models are also dimensionful. For our plots we work 
with units where 6M'^ = 1. 
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Eq. ( 5.15 ), the zero mode satisfies 



Jo 



(0)' 



O AO) , r, .(0) c- AO) 
-JOO/o + "JOl/i + '-'02/2 



A^) I Q flUj , r, A^!) , r, AD) 
Jl + ^OlJo + '^ll/i + '^12/2 

f{0)' c- f(0) I C f{0) I c f(0) 

72 + '^02/0 + '^12/1 + 'i'22/2 



0, 
0, 
0, 



(6.23) 



where the matrix elements Smn are given explicitly by Eq. ( |5.23 ) using the superpotential 
of Eq. ( |6.19| ). It is easy to verify that the normalizable zero mode solution to the above 
system is given by 



\f'. 



(0)N 



(0) 
(0) 



V 



V 

^e-3/2a 





(6.24) 



where A'^ is a normalization constant. This zero mode physically corresponds to changes in 
the size of the extra dimension. 



Example 2 

The second model we would like to analyze has the following superpotential 

W (<I>i, <I>2) = Oi sinh(z^<I)i) + (a<I>2 - h<t'^ 



(6.25) 



where $1 and $2 are again the dilaton and the kink fields. We write the fields in terms 



of their background and perturbations as in Eq. (6.20). As we have explained before the 
parity requirements on the gravitational background a forces us to have an odd parity 
superpotential. To satisfy this, we choose the background kink solution to have odd parity, 
as in the previous example, as well as the background dilaton solution: 

(/>i(0) = 0, </'2(0) = 0. (6.26) 

For this superpotential, the background solutions can be obtained analytically as 

2 , / ai/2y\ 
— arctanh tan , 



h{y) 
h{y) 



tanh(V3a6 2/) . 



The position of the soft- wall singularity in y coordinates can also be analytically determined 
to be 

(6.28) 



Vs 



which is the location in the extra dimension where the dilaton field and warp factor diverge. 

The zero mode solution is given as a solution to = 0. The system of equations 

are given by Eqs. ( |6.23| ). We solve this system numerically using the superpotential in 
Eq. ( 6.25| ) . Parities of the perturbations must be the same as the parities for the background 
fields: 



/r(o) 



/f)(0) = 0, 



/i^o) 



0. 



(6.29) 
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Figure 2: Numerical solutions for the zero mode solutions and for example 2 with 

the parameter values a — 0.5, & = 0.3, a = 1 and = 1.4. The free initial condition for /o(0) is 
chosen to be 1.0. The figure shows that the zero mode solutions diverge in the region close to the 
soft wall, making the solution not normalizable. These solutions have definite parity and we have 
only plotted the z > half. 



where the scaling property of Eq. ( |6.23D allows us to choose f^\^) = 1 without loss of 
generality. If there exists a non-trivial zero mode in this model, then it must possess 
these initial conditions. Therefore, if these initial conditions lead to a non-normalizable 
solution, the only normalizable solution is the trivial one. We find numerically that the 
solution diverges at the location of the soft-wall, as shown in Fig. ^, and the solutions are 
indeed not normalizable. This lack of a zero mode is consistent with our criterion since, 
in this second example, all the scalars have odd parity. Furthermore, the existence of a 
non-normalizable solution to Eq. ( |6.23 ) is consistent with the proof of our criterion, which 
explicitly relied on the fact that perturbations must be small compared to the background 
solution. 

In summary, the model presented here does not contain a zero mode and the physical 



size of the extra dimension is stabilized at the value given by Eq. (6.28). 



7. Conclusions 

Models with warped extra dimensions provide for interesting and rich extensions beyond 
the Standard Model. The classical backgrounds of such models generally contain bulk 
scalars with non-trivial profiles. Stability of the background is an important theoretical 
issue, in particular, a model with a warped, compact extra dimension must have the size 
of the extra dimension stabilized. In this paper we studied perturbative stability of 
scalars coupled to gravity by analyzing the spin-2 and spin-0 spectra. Although our results 
are quite general, we paid particular attention to the case of domain-wall models with a 
soft wall in an AdSs background. These models are interesting to study because, as we 
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have shown, they provide a purely field-theoretic mechanism for compactifying an extra 
dimension. 

We have done part of our analysis using the fake supergravity approach, which has 
been widely exploited to engineer analytically tractable solutions to Einstein's equations. 
This approach has also been used to study stability of models where a single scalar field 
with non-trivial profile is coupled to gravity in a warped background [28, pH, p^]. In this 



paper we have extended these previous studies to the case where an arbitrary number of 
scalar fields couple minimally to gravity. 

The first main result of this work is the derivation of Eq. ( 5.2l| ) , the coupled Schrodinger 



equation governing the spin-0 spectrum of N scalars coupled to gravity. This equation is 
valid for an arbitrary scalar potential, with or without additional brane terms. Following 
this, we specialized to potentials generated using fake supergravity and presented our 
second main result: a system with N scalar fields coupled to gravity and without branes 
has no tachyonic modes in the spin-0 sector. This general result is valid for all types of 
models, where the extra dimension may be infinite in size or finite with a soft-wall, and 
where there may or may not be definite parity. Extensions to models with branes that have 
brane potential terms were also briefly discussed. This result generalizes previous studies 



on stability of models with one scalar field coupled to gravity |28, Using our formalism 



we also studied the case with one scalar field and provided a lower bound on the mass of 



the first spin-0 mode, Eq. ( 6.14 ) 



Our third main result is related to the existence of zero modes for models with 
scalars coupled to gravity. Zero mode solutions in general destabilize the size of the extra 
dimension, and our aim was to determine criteria which guaranteed the absence of such 
zero modes. A general analytic study for models with an arbitrary number of scalar fields 
is rather complex, and we restricted ourselves to scenarios where the extra dimension has 
definite parity. We have proven a criterion that relates the number of zero modes to the 
parities of the scalars and used this result to show that zero modes are absent in models 
where all background scalar profiles have odd parity. We demonstrated this by explicitly 
constructing two domain-wall models with a soft wall, one of which admitted a zero mode 
and the other not. The latter is an example of a model that stabilizes a compact extra 
dimension without using dynamical branes. Whether these models are realistic models 
that can solve the hierarchy problem is a question to be investigated in future work. 
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A. Analyzing the eigenvalues of the Hamiltonian 

Take a system with an arbitrary number of scalars that satisfies Schrodinger's equation, 

H'^ = E'^. (A.l) 
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Below we show the standard result that if one can write H as in supersymmetric quantum 
mechanics |32], 

H ={d, + S){-d, + S)+K^K , (A.2) 

with Hermitian S, the eigenvalues of H are non-negative for vanishing boundary terms. To 
see this we multiply Eq. (A.l) from left with '^'^ and integrate over the extra dimension, 
which gives 

E j dz\^^ = j dz^'fH^ 

= j dz^'t(2 + 5t)(-S + 5)^ + j dz\K^\'^ 

= J dz^\-K + S^){-dl + S)^ + J dz\K^'^ + J dzd,{^\-d, + S)^) 
= j dz\S^\'^ + j dz\K^'^ + j (9^(^+5^) ^ (A.3) 
where we defined 

S = i-d, + S) . (A.4) 

The arrows on the partial derivatives indicate which way they act. Notice that the last 
terms on the last two lines are the boundary terms. Since > and ji^T^'p > 0, for 

an arbitrary ^' 7^ 0, we can immediately see that ii^ > if the boundary terms are zero, 
that is 

^^5^-1 boundary = ^ ^ > , (A.5) 

which is satisfied when either boundary = 

or cS^-l 

boundary — 0. In fact the requirement 
that a Hamiltonian is a self-adjoint operator for a physical problem already forces these 



boundary terms to vanish. We explicitly verify this for our particular models in Section 5.2 
Note that these conclusions also hold when K = 0. 
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